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r-| ■ Abstract 

p • We study the energy and nuclear A dependence of the hadronic production of 

^ ■ heavy quarkonia. We review theoretical ideas which have been put forward, 

seeking a consistent global picture reconciling the large effects in quarkonia 
^ ! with the small nuclear effects observed in continuum Drell Yan production. 

The data indicates that shadowing or leading twist modifications of parton 
distributions can be ruled out as explanations, leaving higher twist energy 
loss. From general principles the maximum allowed energy loss of partons 
traversing the nuclear medium can be related to the parton transverse mo- 
menta. We then show that the experimental data on nuclear suppression of 
charm- and bottom- onium for large Xi;' is consistent with this effect: using 
the observed transverse momenta to bound the xp dependence in an almost 
model independent manner generates a relation that practically reproduces 
the data. Several prediction are discussed; the dependence on xp as 1, 
and large and small cuts, can be used to discriminate between quark and 
gluon induced effects. 
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1. INTRODUCTION 

Propagation of quarks, gluons and hadrons through nuclear matter is cur- 
rently a subject of intense interest. It is expected that the study will teach us 
much about the interplay between perturbative and non- perturbative QCD. 
An important experimental discovery by the Fermilab E772 experiment [1] 
is the suppression of charmonium and bottomonium production in pA colli- 
sions in comparison to the production rate in pp collisions. A similar effect 
had also been previously reported by Badier et al [2] and Katsanevas et al. 
[3]. The data strongly contradicts a widespread theoretical expectation that 
at high energies the nuclear medium should have negligible effect on heavy 
quarkonium production. The observed suppression has direct implications 
for the use of onium production as a signal for quark-gluon plasma forma- 
tion in heavy ion collisions. It has generated much controversy [4-7], and 
raised the possibility that the observations represent a serious challenge to 
theory. 

A common theoretical prejudice that suggests negligible nuclear suppres- 
sion is the following. One can argue that the characteristic separation of the 
charmed quark anti-quark pair produced by a partonic interaction is of the 
order of l/rric, which is equal to (1.5 GeV)~^. Invoking a geometrical cross 
section in the spirit of color transparency, the attenuation cross section of the 
charm pair might be as small as the order of (1.5 GeV)^^ , which is about 0.2 
mb. With this cross section, a typical survival factor of 0.97 is obtained for 
a large nucleus of diameter 10 fm. This is a very small effect in comparison 
to 30 to 60 % suppression seen in the data. 

Of course, such a method of estimating the nuclear effects applies (at 
best) to the propagation of a well-locahzed, relativistic, color-singlet charm 
pair. But invoking color transparency for the actual production is probably 
not correct, since the kinematics of the events are highly inelastic, and lack 
the usual conditions of exclusivity that color transparency arguments should 
assume. Color transparency is not expected to occur if the coherence of a 
system is broken, for example in the case when uncontrolled inelastic color 
flows are summed over in a semi-inclusive production. Moreover, in the initial 
state the gluons can lose energy by interacting with the nuclear medium, and 
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there also is the hkehhood that the charmed pair is produced temporarily 
in the color octet (rather than singlet) state. This does not mean that the 
suppression is of the expected size: indeed, the huge magnitude of the effect 
has been quite mysterious. 

Besides the outstanding puzzle of quarkonium suppression, the general 
problem of parton propagation in the nucleus is becoming rather important. 
There is great interest in understanding the role of quantum mechanical 
coherence of QCD interactions in the nuclear medium. Depending on the 
experimental circumstances, the same arguments leading to coherent sup- 
pression of energy loss in high energy QED interactions - the so called LPM 
[8-10] effect - can be applied, leading to many interesting predictions. In the 
LPM effect, a concept of "formation time" Tform of quanta during a high en- 
ergy interaction is considered. The order of magnitude of the formation time 
for a quantum of mass m, carrying energy E, and transverse momentum kx 
with respect to its progenitor is 

Tform ~ E/{k^ + m^) (1) 

The LPM effect in QED suppresses bremsstrahlung due to destructive in- 
terference between emissions occuring over an interaction time Tint < Tform- 
A recent paper by Gyulassy and Wang [11] studies the effects of multiple 
scattering in perturbativc QCD. Although the details of coherence are con- 
siderably more complicated due to the non-Abelian color algebra, the basic 
features of the LPM effect are not changed by this study. 

The concept of the formation time leads to a complementary concept 
of formation rate Fform- The formation rate is simply the inverse of the 
formation time, 

Tform ~ (4 + m'')/E (2) 

Examined in this way, there is an interesting possibility that the coherent 
parton formation rate could be "tuned" by selecting signals with various 
masses, energies and transverse momenta. Is it possible, for example, to 
increase particle formation by selecting events with larger k^l As we will 
show below, the trends in the data indicate possible observation of a dramatic 
"anti-LPM" enhancement of parton emission due to increased formation rate 
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associated with large k^. The possibihty of tuning the formation rate leads 
naturally to a number of interesting signals which can be experimentally 
tested. 

2. REVIEW 

In view of these points (and the historical difficulty of understanding 
charm production in general), we will approach the problem of onium sup- 
pression in a methodical manner, and try to limit our assumptions carefully: 

i) Electroproduction experiments have convincingly shown that the process 
of parton fragmentation is negligibly affected by nuclei. The E772 experi- 
ment itself also found that continuum Drell-Yan muon pair production shows 
no significant nuclear effect. The data for continuum dimuon A~^da/dQ'^, 
for example, shows little nuclear effect in the Q'^ regions both above and 
below the suppressed onium regions. These observations are consistent with 
conventional factorization. They seem to significantly limit the amount of en- 
ergy loss one can attribute phenomenologically to parton propagation inside 
the nuclear medium. 

However, we note that both the bulk of fragmentation and the Drell 
Yan continuum are mainly probes of the quark- antiquark channel. They say 
little about production via gluons. The nature of gluon and heavy quark 
propagation in nuclei is not clear. 

ii) A very important clue to the charmonium data is given by the fact that 
the data's X2 dependence does not scale. Even if one creates an ad-hoc gluon 
or quark-antiquark distribution for one experiment on charm, it does not 
reproduce the data for bottom. Furthermore the data docs not scale with 
X2 as wc change the incident beam momentum from 200 to 800 GeV [12]. 
This eliminates the possibility of gluon shadowing or more general EMC-type 
effects on the parton distributions as a dominant mechanism. Quite recently, 
Banesh, Qiu and Vary [7] (BQV) claim that the onium production at large 
xp is dominated by quark annihilation, and that shadowing in the small- 
X2 quark distributions accounts for the dominant part of the suppression. 
Of course, shadowing makes a contribution, but the magnitude of small x 
shadowing observed in deeply inelastic scattering is not nearly large enough 
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to account for the onium suppression. Postponing a detailed discussion to 
Section 5, we do not believe that the BQV analysis actually support the 
claim. The same objection based on factorization must be raised here: the 
onium production data shows that factorization does not hold, ruling out 
this mechanism as a model for the full effect. 

iii) Since factorization is a leading twist prediction for this kind of semi- 
inclusive production, and the same sort of experiment has confirmed its use 
and understanding elsewhere, the suppression effect must be a higher twist 
one. The bizarre thing is that it so drastically affects the onium production. 

iv) The quarkonium suppression seen in the E772 data can be questioned 
as to whether it might be an instrumental effect. In fact, certain regions 
of the data's transverse momentum spectra have not been reported due to 
questions about the acceptance [13]. We re-examined previous data of Badier 
et al [2] on J/il^ production on nuclear targets in studying this. We find that 
rather than contradicting the E772 experiment, the Badier et al data seems 
to confirm the trend. A signal for suppression of onium is real and more than 
ten years old. 

We now summarize some theoretical ideas, which we consider as "spare 
parts" that might be assembled in a new order to understand the puzzle: 

a) It can be argued that heavy quarkonium, as a non-relativistic bound state, 
might have a large nuclear matter interaction cross section that is set by the 
full onium size and binding energy rather than the mass. The J/%1) bound 
state diameter, for example, is 3 to 5 times the charmed quark compton 
wavelength, leading to an estimated cross section as large as 2 to 5 mb. A 
major problem with this proposal is that the time scale for formation of a 
bound state onium is very long compared to the time scale for crossing the 
nucleus. This is aggravated by the fact that the Lorentz boost stretches any 
effective time scale enormously, so that a J/ip does not become a J/ip until 
it is more than 100 Fm away! We do not pursue this idea further. 

b) There remains a realistic possibility that interaction of incoming and out- 
going colored partons could cause them to lose energy. Gavin and Milana 
(GM) [5] observed that even a small shift in the Xp value of the incoming 
partons, assumed by them to be gluons, could lead to a numerically large 
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change in da/dxp because of the rapid variation with x of the gluon distri- 
bution functions. 

For example, supposing the gluon distribution to go like (1 — x)^, then a 
gluon taken from the gluon distribution aX x + 5x is less likely to contribute 
by the ratio {1 — x — 5x)^ / (1—x)^ ^ 1 — 55a;/ {1—x). Even if 5x is small, this 
creates a strong kinematic suppression as x approaches unity. For reference 
we will call this the GM mechanism. This effect was assumed to occur for 
both initial and final state propagation, with onium production dominated 
by color octet components. 

To get a large enough shift in the x value, GM also proposed a rule for 
energy loss which is of the "higher twist" type; their proposal for the energy 
loss is ^E/ E = cx\A}l^ IQ^ ^ where c is a color dependent constant that could 
be adjusted to fit to the data . At first the higher twist character of the GM 
rule seemed to put it into a phenomenological limbo of the uncalculable, a 
thing which could be neither verified nor disproven with current theoretical 
knowledge. 

c) This approach was countered by Brodsky and Hoyer (BH) [6] who argued 
that the dependence on energy of the GM formula violates general principles. 
BH went on to claim that there exists an upper bound on the energy loss 
for a parton propagating through nuclei. This upper bound, in the spirit of 
LPM, is obtained from rather general considerations; the BH rule is Ax < 
A;^L^/2£', where is the transverse momentum change in the collisions 
giving the energy loss and La is the target length. This rule also is higher 
twist, but the uncalculable higher twist mechanism of Gavin and Milana has 
a limit that contradicts this relation, allowing A£' to go like E at fixed Q^. 
The contradiction between the two formulas is numerically important. Using 
their own bound for the energy loss and a value for k^— 300 MeV, Brodsky 
and Hoyer dismissed the Gavin and Milana proposal as insufficient to explain 
the data. 
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3. ASSEMBLING THEORETICAL SPARE PARTS 

Considering these ideas, let us first observe that the GM mechanism is 
(a) quite reasonable and (b) logically independent of the model used for the 
energy loss. A first task, then, is to determine whether some energy loss, 
however it occurs, can explain the data. This has already been answered by 
GM who were able to fit the data. There is no reason, then, to rule out gluon 
energy loss as a mechanism, but we agree that one should have a consistent 
framework to represent it. 

A second observation is that the BH expression for energy loss can be 
viewed a testable hypothesis, that is, as a model. The rule is claimed to be 
an upper bound, so any observed energy losses must lie inside an envelope of 
values it specifies. This prompts us to examine the experimental data, using 
the GM hypothesis, and test for the general mechanism of energy loss by 
seeing whether or not the data obeys the claimed bound. Our goal here is to 
check the energy loss proposal without getting snarled into model dependence 
of the energy loss formula. 

Third, given the fact that Nature tends to dissipate energy rather maxi- 
mally, we can try saturating the bound and examining whether we obtain a 
prediction close to the behavior shown in the data. Actually, by examining 
the bound more closely, we find it needs modification from mass effects and 
also by an unknown dimensionless factor. We consider our additions to be 
modest corrections. The trends in the data are not strongly dependent on 
the unknown prefactor. 

The arguments leading to the BH bound, its modifications, and its appli- 
cation along with the GM mechanism are reviewed in the next section. We 
find that we can test for energy loss by examining the detailed xp distribu- 
tion using the /cy distribution as an input. This is a much more powerful test 
than simply looking at overall production rates. The procedure works quite 
well in the upsilon meson case where the kx distribution has been measured. 
We find that the suppression seen in the xp distribution is well within the 
limits imposed by the bound on energy loss. Even more importantly, the 
dependence of the data on Xp and A actually tend to parallel the bound. We 
consider the agreement of a general bound and a hitherto unnoticed pattern 
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in the data itself to be practically model independent evidence that the basic 
culprit in nuclear quarkonium suppression is gluon energy loss. 

While the bottom quark case checks well, a definitive application of this 
result to E772 charmonium production is not yet possible because of prob- 
lems due to experimental acceptance [13]. The transverse momentum distri- 
bution for the charmed case is known only for kx less than about 2 GeV. We 
can nevertheless apply our formalism in reverse and use the experimentally 
observed xp distribution to put a lower limit on the distribution of char- 
monium. This prediction can be tested in future experiments. Our results 
show that the lower limit on the /ct distribution is roughly the same as the 
corresponding observed distribution for the case of bottomonium. 

3.1 THE ROLE OF GLUONS 

Although we will present evidence that the GM mechanism is at work, this 
connection seems, at first sight, not specific to the production mechanism. 
There is a problem due to the generality of the BH bound. The BH bound 
assumes too little; one might say it tells us that we are seeing the uncertainty 
principle at work, a fact of limited usefulness. One knows energy was lost 
but does not identify the underlying subprocess with this mechanism alone. 

The surprise of the onium suppression is contained in the large magnitude 
of the effect. One of our main points is that the quarkonium suppression 
(which has gotten so much attention) is directly related to the transverse 
momentum in the data (which has gotten very little attention). The average 
transverse momentum squared < k'^ >a is rather large;[] for Tungsten the 
< >A for bottomonium is greater than 2 GeV^ rather than 0.1 GeV^. As 
a function of A, the integrated Tis data is well described by [1,14] 



The same experiment found the A dependent part of the Drell Yan continuum 



similar observation on the largeness of nuclear interaction induced tranverse mo- 
mentum in dijet production has been made by T. Fields and M. Corcoran (to appear in 
Proceedings of EPS Conference^ Marseille 1993). We thank Tom Fields for informing us 
of this. 




(3) 



< k"^ >A to be about 10 times smaller. As we will show, much can be 
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predicted simply knowing the < >a- 

What, then, is causing the rapid variation with A of the quarkonium 
< k'^ >A^ Variation of mean transverse momentum with A has also been 
seen in other experiments studying nuclear dependence of dijet and dihadron 
production [15,16]. Physically this effect can arise from multiple scattering 
in the nuclear medium [17]. The data's A^^^ dependence indicates a target 
length effect, consistent with the energy loss mechanism. The same data 
indicates that the charged partons contributing to Drell Yan are acting dif- 
ferent from the partons making the heavy quarkonia. Drell Yan precludes 
any strong effects for the quarks, then! This suggests that gluons, which 
give dominant contribution to heavy quarkonia production, are scattering 
more in the transverse direction and are losing more energy than quarks. An 
alternate possibility, which is not strictly ruled out by the current data, is 
that the final state interactions of the heavy quarks are responsible for the 
large nuclear effects. However, experience with QED bremsstrahlung and 
perturbative QCD, where a quark mass acts to substantially cut-off mass- 
less vector emissions, makes it hard to believe that heavy quarks could lose 
energy so much faster than light quarks! In Section 5 we show that analysis 
of data at large xp will also be able to discriminate clearly between these 
two possibilities. 

We believe that further data and not theoretical arguments will play the 
most important role in finally pinning down all the unknowns. 

For better or worse, the dynamics of gluon channels producing onium 
have never been clear, and a number of issues, including "intrinsic charm" , 
have been created to address the problem. One cannot assign a perturbative 
overall normalization to the gluon channels safely. It is not known, e.g. what 
fraction of the time a color octet qq state is produced rather than a color 
singlet. For this reason, we have arranged our calculations so as not to base 
them on normalization factors. We will postpone to Section 5 a discussion 
of the interplay of subdominant quark channels with the gluon channels. 

3.2 ENERGY LOSS RATES 

We first review and expand on the argument of Ref. [6] to obtain an 
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expression for an upper bound on energy loss. Assume a parton propagating 
through the nuclear medium in the +z direction with energy E — xiEp, 
where Ep denotes the energy of the incoming proton. This parton loses an 
energy XgE by emitting a gluon in the presence of a source as shown in Fig. 
1. This, and any secondary scattering, must occur in the volume of a nucleus. 
The finite size of the nucleus introduces a distance scale Ax, which will be 
used with the uncertainty principle to bound the energy loss. 

Components of momentum q exchanged with the source can be related 
to the source spatial size Ax by AxAq > 1. The "— " component of the 
exchanged momentum is easily probed. Letting the four momentum squared 
of the final state be given by M^, then = + 2p ■ q. We solve for 
q- fti Aq- fti MV(2p+) - q'^/{2p+), where p+ = {E + p^)/y/2 

Assume for definiteness now that we have two identical particles of mass 
m in the final state. Their momenta in (+, T, — ) notation can be listed as 



Pi = 
P2 = 



XgE, Ut, k^/{2xgE) 
'{l-Xg)E, -It, kl/{2{l-Xg)E) 



where /c^ = + is the "transverse mass" . The final state mass is 



X n 



+ 



X n 



+ 2 



X{] 



Xfi 



(4) 



(5) 



Assuming k^ « M^, then the kinematics forces either Xg « 1 or (1 — 
Xg) « 1. For gluons emitting gluons the two situations are physically 
identical. Thus it is a good approximation to take 



k% + 



Xn 



(6) 



Combining the above with Aq — (M^ — q^)/{2p'^), we have 



1/Aq- ^ 2p+Xg/(4 + m^- Xgq^) , (7) 

This can be combined with the null plane uncertainty principle Ax'^Aq~ > 1 
to give 
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2p+Xg/{kl + - Xgq^) < Ax+ , (8) 

The meaning of Ax"*" is the hght-cone "time" within which the events 
occur. This is clarified by a space-time cartoon (Fig 2). From the figure, 
due to the finite size of the nucleus, the Ax"*" for events causally propagat- 
ing along the future light cone obeys Aa;+ < \/2La, where La is the rest 
frame length of the nucleus. Note that non-relativistic propagation would 
allow Aa;+ to become indefinitely large if the parton stops inside the nucleus. 
This possibility is irrelevant to the discussion, although it would weaken the 
bound. 

BH did not include the effects of and q^. Ordinarily such terms are 
negligible in comparison with large momentum transfers. However in this 
case we see that they are not necessarily negligible in comparison with k"^. 
Setting and to zero in (8), we obtain the result of BH 

Xg < LAkl/{2E) (9) 

using p"*" = s/^E. The energy loss is then given by AiJparton ~ XgE — XgXiEp, 
leading to Axi — XgXi < LAk^/2Ep, as found by BH. 

Restoring the dependence on and q'^, we find: 

^ {LA/2E){k'^ + m') 

l-\q'\LA/i2E) ■ ^'"^ 

At this point, if we know \q^\LA/2E « 1, then we have a BH bound with 
only the modification of the "transverse mass''. For our further discussion 
we will assume this limit. 

3.3 BACKUP 

Are there any important loopholes in the bound? We note the following: 

i) The neglect of a final state mass for the parton is dangerous. In the BH 
calculation, using kx values around the traditional values of 300 MeV, then 
k"^ = 0.1 GcV^, which is much smaller than any value one would use for 
an effective m^, even considering "massless partons". The bound is quite 
sensitive to this. We will use k^ values obtained from data. This is an 
important detail for application of the bound, but not for its concept. 
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ii) The bound, while invoking general kinematics and quantum principles, 
nevertheless depends on the dynamics which assumed the production of two 
particles per vertex in the final state. It is similar in spirit to electrodynam- 
ics, where the LPM effect [8-10] provides the classic prototype for coherent 
suppression of energy loss in a finite density target. 

To see this, we present a "back of the envelop" discussion of the physics 
of LPM, abstracted from Feinberg and Pomeranchuk [10]. The kinematics 
assume nearly forward scattering of an electron (mass = mg, energy =Ee) 
in a classical medium, with emission of a bremmstrahlung photon carrying 
energy uj. All transverse momentum are assumed to be very small - this is 
an important point. For forward scattering, one finds a spatial momentum 
transfer g|| given by 

q^^^^El-ml-^{E,-u;y-ml-u; (11) 

so that 

Agii ^ mlu:/{2El) (12) 

This momentum transfer determines a coordinate space region of longitudinal 
length r ~ 1/Ag|| ~ 2El/{rn1uj). Suppose that the length is so large that 
the electron undergoes several multiple scatterings. From multiple scattering 
theory, the electron does a random walk with mean square scattering angle 

where Eg is an energy scale and L a radiation length. Note that 9l is propor- 
tional to the distance travelled r. The LPM argument, which is semiclassical, 
observes that coherence can be retained in bremmstrahhmg emission over a 
cone given by an angle set by the boost parameter 7 = E/m. Most photons 
are emitted with angle < I/7 relative to the moving source. Coherence 
over the emission becomes crucial if the multiple scattering angle Og is bigger 
than 9^. This condition is 

¥\'-r>-m\ (14) 



,EeJ L \E^ 

which upon inserting r ~ 2El/m1uj becomes 



^ ^ n4z (i«) 
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Emissions satisfying this bound are suppressed by destructive interference. 

The emphasis in the LPM analysis is on a region where the photon forma- 
tion rate is small compared to the collision rate. But then it follows that the 
process is exquisitely sensitive to the scales setting the rate Fform = k\j2uj. 
By adjusting k\ we can evidently scan across a range of formation times, and 
turn the LPM suppression on or off. We will apply this observation to the 
quarkonium suppression in the next section. 

Continuing, the dynamics of QCD has features which are different from 
QED. In QCD an incoming gluon can split into three gluons at a single in- 
teraction due to the perturbative 4-gluon vertex (Fig. (3)). This upsets the 
QED-based argument. We have worked through the kinematics of three par- 
ticles in the final state, verifying that some regions reproduce the formation 
time argument while other regions exist which do not. 

For definiteness consider Fig. (3), in which an incoming gluon splits into 
three gluons carrying momentum fractions xi, X2 and x-i — \ — x\—X2- The 
momenta are given by 

p. = Up, kr,, |i) 



X2P, kT,2, 



X3P, kT,3, 



Jvj-< 



4^ 



2X2PJ 



krp 



2xsPj 

in T, — ) notation, assuming and — > 0, and using massless gluons. 
The q" momentum is simply given by the sum of the p~: 

1.2 l2 7,2 

fvrpi tvrp2 I^T,3 

q = _ ' + _ ' + 



2xiP 2X2P 2X3P 
Applying q~ > 1/L, where L is some interaction length, we can bound this 



sum. 



The "formation time" is Tform — ^/q ■ We see that its inverse is the sum 
of three inverse formation times, 

l/Tform = 1/ri + I/T2 + I/T3 (16) 
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where each formation time (up to trivial factors) equals the usual definition 
l/Ti — Ti 2xiP/k^^. The interpretation of (16) is of course simpler in terms 
of rates: the total formation rate Ttot is the sum of the three individual 
formation rates. 

From the uncertainty principle the formation of the final state requires 

^totL ^ 1 

or 

La>^ (17) 

The finite value of La means that not all the formation rates can be too 
small or else they will destructively interfere. 

Unlike the case of 2-body formation - where the creation of one particle 
implies the other - when we have a 4-point interaction, there is more than 
one independent formation rate. In the relation (17) the biggest formation 
rate wins. That is, in the region 

h.2 7,2 7,2 

2xiP 2x2P 2X3P 

then the creation of parton "1" dominates the issue of coherent formation. 
The formation of partons "2" and "3", while occuring less rapidly, is trig- 
gered by the formation of parton "1" at the 4-point vertex, but they are not 
separately resolved at this point. 

Physically, here is what happens. An incoming gluon can be disrupted by 
the source into a small x gluon (say) , while two larger-x comoving gluons are 
simultaneously created. The uncertainty principle applies to the smallest-x 
gluon which cannot be produced too slowly. But understanding the coherence 
and quantum mechanical resolution of the smallest-x gluon says nothing about 
resolving two other gluons into separate components. A subsequent hard 
collision (or similar independent time scale) is needed to resolve them. 

Naturally we have a probe of a gluon's momentum when a heavy quark is 
produced. Suppose, between X2 and x^, we detect a quark carrying X2 ~ Xp] 
what fixes the value of x^l It is fixed by detailed dynamics, not general 
principles. This situation is unprecedented; it indicates the possibility of 
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energy losses not following the rules of QED. We conclude that the usual 
QED-LPM arguments are inadequate for a quantitative analysis of energy 
loss in QCD. 

This is a real loophole; how big are its effects? The four gluon vertex 
is higher order in perturbation theory, being of order g'^, but cannot be 
negligible because it is necessary for the gauge invariance of the theory.0 
Moreover, the problematic integration region is comparable or larger in size 
to the region being discussed in the QED formation time arguments. Yet our 
analysis of the dimensionless 4-gluon vertex emission does not introduce any 
new "large" scales beyond L^, and 1/ E already present. By dimensional 
analysis, then, something like the LPM analysis should survive after detailed 
integrations and combinatories over dimensionless quantities are evaluated. 
We will accept it for this study. We believe that further work could show 
that the bound might be multiplied by a dimensionless factor we estimate to 
be a few units. 



4. APPLICATIONS 

We now turn to applying the energy loss bounds as practical tools. The av- 
erage value of k"^ in our formula represents the transverse momentum caused 
by scattering. We assume that this is the difference of intrinsic transverse 
momentum < k'^ >int, of the order of 0.5 - 1.0 GeV^, and the observed value. 
A more crisp definition can be given for "leading twist" reactions but does 
not exist for power suppressed processes. To proceed orderly we separate 
the initial state energy losses of gluons from final state ones of heavy quarks. 
We present bounds based on initial state gluon energy losses only. We will 
show that we produce a trend that is strikingly consistent with the data. We 
discuss separating initial from final state effects in Section 5. 

First we examine an analytic estimate of the energy loss. We assume 
gluons are bremmed off with a distribution 

dN 

^It has always been worrisome that the 4-gluon vertex has produced few qualitative 
effects in high energy phenomenology. 



15 



where Xmax is given by the bound (9) . Since we are in quasi- nonperturbative 
region the details of fg are unknown. So long as f{xg) is peaked at small Xg 
but regular at Xg — the details turn out not to matter much. 

The effective gluon distribution G{x) due to the shift from energy loss is 
given by 

G{X) = / dXgf{Xg){l -X- Xgf 

JO 

Using for example f{xg) — 1/xg, x^in < Xg < Xmax- then we are interested in 
the limit x^ax << 1- A series expansion can be obtained by integrating by 
parts 



^■^111 ax 

dXn 



f ... 



Hj<uU fj 



X — Xg)^ = 111X^(1 — X — Xg) 



max 
^min 



+ 5 / dXg\n(Xg){l — X — Xg) 



The second term can be integrated by parts again leading to an asymptotic 
series of any order desired. The first term is approximately 

5xti 



G ~ In (Xmax/a^min) {I - x)- 



^ — max ^ 
1 — X 



to first order in Xmax and dropping terms proportional to Xmin- The logarithm 
is slowly varying in both Xmax and the infrared cutoff Xmin and will be ignored. 
One sees the GM mechanism emerging in the power series expansion: the 
effects of a small x^ax get big as x — > 1. 

Now suppose that G is used in a calculation of cross section, namely 

XiG{xi) X2G{X2) 



dk^dxidx2 dk^dxidx2 
where da is evaluated at the parton level. Setting up the kx integrals we 
have 



, = / d^kTXiG{xi) X2G{x2) ( 1 - ^^r-^A 

dxidx2 J dkj'dxidx2 \ 2E[1 — xi] 



The second terms contains the nuclear effects; we have inserted the bound 
Xg = k'^LA/2E. Since the integrand is proportional to k^, we can do the 
integral to estimate the A-dependent correction as 

1 daA _ dai / 5 < k^ >a La ^ \ 
Adxidx2 dxidx2 \ 2E{1 — xi) I 
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The correction has a size set by < >a- For nuclei La ~ 1-2 Fm A^^^, and 
the data for Ti^ production gives < >a~ 0.16^^/^ GeV^+ <kT>D- The 
kinematics of producing quarkonium with invariant mass and momentum 
fraction xp requires xp — X1—X2, Q^/s — X1X2, which in the hmit Q'^/s « 1 
gives Xi? ~ Xi- Then for the estimate for T15 production 

1 daA/dxp , 5 < A;^ >a La 
Adai/dxp 2E{1—xf) 

This crude estimate does surprisingly well. Take for example A = 184 for 
tungsten, and xp = 0.5. Then the effects of including energy loss lead to 
about 60% of the events compared to a calculation neglecting energy loss. 

To check the analytic estimate wc did some numerical calculations. We 
wish to compare the experimental trends in the data's x-dependence with its 
kr dependence. We will use the experimental value of < k'^ >a — < k"^ >int, 
where < k!^ >int is the intrinsic value, to calculate A,ti. This can then 
be used to calculate numerically the shift in the gluon distribution functions 
due to energy loss, thereby yielding the x dependence in the nuclear medium. 
Thus we take from the data 

Axi < La{< kl>A-< kl >int)/{2Ep) . (18) 

where we set < k'^ >int to be equal to 0.91 GeV^ in analogy to Drell Yan 
[18]. 

The x-integrated kx distribution of the quarkonia has been parametrized 
experimentally [1] by, 

1 6 



1 



1 + {kr/poy 



(19) 



where ^ {A) is an A dependent normalization factor and po is also A depen- 
dent. The average value of p^, defined by. 



^2 ^ _ lo dkj.kj.fA{k^) 



<^T>A- J^dklfAikl) ' 

is equal to Po/^- As given in Ref. [1], the values of po foi' are 

Po = 2.78 for Drell Yan , 
Po = 3.22 for T . 
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As discussed in the introduction, the values of po are unfortunately not 
available for charmonium. The details of the dependence are not our 
object here, and we would prefer to insert them from data. To proceed 
with charmonium we will assume that it has the same transverse momentum 
dependence as bottomonium-an assumption which can be relaxed when data 
is obtained. The value of < >a for the case of bottomonium increases as 
0.16A^/^ [14]. The experimental fit to this data for the case of T is given in 
eq. (3). 

We next need the production cross section in terms of the parton dis- 
tributions. We assume that gluons give the most important contribution to 
quarkonium production [19] for moderately large Xf - this is discussed in 
detail later. The cross section integrated over transverse momentum is then 
given by 

da xiX2 ^, ^ ,a39^^%Q^) 

- ' -G{xi + Axi)G{x2) 7^^^- (21) 



dQ'^dxp X1 + X2 

We apply the same procedure to the Drell Yan continuum dilepton data, 
substituting quark and anti-quark distributions for the gluon distributions. 
The transverse momentum is quite different, but has been re-fit to match 
the data. The overall normalization is not relevant because we report ratios 
of nuclear targets to the proton. This procedure, repeated for each nucleus, 
gives us definite predictions for the x- dependence of the experimental data, 
which is discussed in the next section. 

5. RESULTS 

In this section wc discuss the results of our simple GM energy loss cal- 
culation combined with the modified BH rule. The xp dependence of the 
ratio [do dxp)j doi)! dxp) extracted by using the experimentally measured 
/ct dependence is shown in Figs. (4-7) for the cases of Drell Yan, Ti^ and 
charmonium respectively. For the bound wc set Axi (Eq. 18) to its max- 
imum allowed value. Wc note that the experimentally measured points are 
well within this theoretical limit; as mentioned earlier, the trend in the data 
is to run parallel to the bound. For the case of charmonium we have taken 
the transverse momentum distribution to be the same as for bottomonium. 
We see that this choice fits the charmonium xp dependence and therefore 
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gives the minimum value of the transverse momentum for the case of char- 
monium. The numerical results also include shadowing besides energy loss. 
We assumed the following functional form for the ratio it!shadowing of nuclear 
to deuteron quark distributions due to shadowing, 

i^shadowing = 0.809 + 0.261 exp (-^2 - 0.00526^1^/72^2) 

which fits the structure function data [20] with ^^/(degree of freedom) = 
0.86. The shadowing for the case of gluons is included by assuming that it is 
the same as for quarks. The dashed lines in fig. (4), (5) and (6) represent the 
results without including energy loss. We see that although the Drell-Yan 
case is well reproduced by shadowing, charmonium and bottomonium data 
cannot be explained purely be shadowing, but require the addition of gluon 
energy loss. 

Our results lead to several experimentally checkable predictions: 
Small kt 

First, the conventional A"' analysis can be examined bin by bin in kr and 
xp. Generally speaking, the energy loss picture is distinguished by producing 
the largest suppression in the largest /cy regions. This does not mean that 
small kr is totally safe, because rescattering can feed particles back in to 
this region. However we expect this to be controllable and therefore predict 
suppression increasing monotonically with kx at fixed xp ot integrated over 
xp- This implies that the suppression should be reduced if we consider the 
bin with kr less than about 1 GeV^, where the dominant contribution comes 
purely from the primordial tranverse momentum. Our estimate of the xp 
distribution for different kx cutoffs is given in Fig. (8). The perturbative 
part of transverse momentum is again calculated by subtracting the intrinsic 
contribution from the observed transverse momentum. At low transverse 
momentum, < 2 < >int we calculated the limiting value of Axi by 
setting k^ equal to the intrinsic value of 0.91 GeV^. 

The resulting curves, in Figure (8), show the ratios of the cross section as 
a function oixp for various values of kr- The curves are on a log plot, because 
the overall normallization is not being predicted. A shift in normalization N, 
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translating the plots up or down, can be considered a free parameter. Our 
object is the shape of the curves, which clearly evolves with the kr cut. For 
each region the experimental data should lie above the corresponding curve 
and follow the trend indicated in Fig. (8). 

Large xp 

The idea of dominance of gluons in heavy quarkonia production is hardly 
new but has been raised again recently by BQV [7]. They consider the lowest 
order pcrturbative subprocess cross sections for producing a quark- antiquark 
pair with invariant mass Q^: 
quark- antiquark channel: 

gluon channel: 

where 7 = Aml/Q"^. Convoluting these cross sections with standard parton 
distributions, BQV claim that the quark initiated process dominates over 
the glue-glue one for xp ~ 0.5. This conclusion is based on the perturbative 
normalizations given above, and the fact that the quark distributions fall less 
rapidly with x than the gluons. Quarks are forced kinematically to dominate 
as xp goes to 1. 

We agree with this in principle, but disagree that the crossover point 
can be given by the Born term calculation. Long experience with detailed 
calculations of quarkonium production at high energies favors gluons over 
quarks, indicating phenomcnologically that the pcrturbative normalizations 
arc not to be trusted too literally. As already noted, the data does not allow 
the option of ascribing the quarkonium suppression to quark channels while 
simultaneously accommodating the dilepton continuum. Moreover, the Born 
term docs not even give the Drcll Yan cross section correctly; for a long 
time it has been known that a "K-factor" of about 2 is needed to fit the 
data. Current understanding of K-factors is that they summarize higher 
order corrections from initial and or final state interactions. One cannot 
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assume the K-factors cancel out in nuclear ratios: what is relevant is the 
relative amount of quark and gluon contribution in each target. 

There are several ways to proceed. One can estimate the crossover be- 
tween quark versus gluon dominated quarkonium production with the K- 
factor method. Using Ref. [21], in the initial state interaction between two 
gluons we find a K-factor which is bigger than the annihilation K-factor by 
Nc/Cf — 9/4. With this (crude) estimate, the crossover point for quark an- 
nihilation channels over gluon channels at 800 GeV can be estimated to be 
0.65. One may also treat this as an adjustable parameter to be determined 
experimentally once more data becomes available. 

The value of large xp is that one can experimentally "tune" the produc- 
tion process to favor quark initiated reactions. We have already noted that 
the Drell-Yan data indicates almost negligible energy loss, and smaller trans- 
verse momentum, for light quarks compared to gluons. If this is correct, then 
as a: is increased above the crossover point the suppression in the nuclear 
medium should diminish. In Fig (9) we present a calculation illustrating 
the effect. This calculation was performed by using the Eichten et al. [22] 
parametrization of the quark and gluon distributions. If the final state ef- 
fects are neghgible than the data should show a sudden change in the current 
trend and an enhancement at large xp- This is a dramatic signal meriting 
a careful search. For this calculation, we used the Born term cross sections 
modified by a relative normalization of 9/4 for the gluons to the quarks. 

Our approach has combined theory with empirical patterns taken from 
the data. One could ask why light quarks do not come close to saturating 
the energy loss bound while apparently gluons do. The answer is, we don't 
pretend to know. In the same vein, one can ask whether final state heavy 
quark energy losses should have been included. The answer is, the data 
does not indicate that significant energy loss from the heavy quarks needs 
to be introduced. Nevertheless, toward developing a truly model indpendent 
procedure, let us observe that the limit xp — > 1 plays a key role. Suppose 
the up-turn as xp — > 1 predicted above does not occur, even at such large 
values oi Xp that we know the production is quark dominated. Then the 
suppression of heavy quarkonium compared to Drell-Yan production must be 
due to the heavy quark interactions above. Similarly, comparing experiments 
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with different beams - and especially pion beams known to be richer in 
quarks as x — > 1 - has the potential to help separate the final state from 
the initial state effects. The Badier et al. data on pion initiated reactions [2] 
is consistent with this trend. We recommend using the data itself to separate 
the issues in a systematic way. 

Certainly the effects discussed here have a direct bearing on the use of 
charmonium or quarkonium suppression as a probe of quark- gluon plasma 
formation at RHIC. Certainly a more thorough theoretical understanding 
and further experimental investigation of the phenomena is required before 
firm conclusions could be drawn from quarkonium production in heavy ion 
collisions. 

Note Added: After this work was completed we became aware of a recent 
paper by M. S. Kowitt et al., Phys. Rev. Lett. 72, 1318 (1994), which 
has extended the experimental data on Xi;' to larger values then was avail- 
able previously. Their results show that the ratio of nuclear to Deuterium 
production starts to rise up considerably beyond xp — 0.65 but then falls 
again around xp — 0.9. Except for the point at xp — 0.9, this data seems 
to support our picture. Comparing the data to our fig. (9), this suggests 
the idea that both light and heavy quarks lose negligible energy compared 
to gluons may be correct. 

Acknowledgements: This work was supported in part under Department 
of Energy Grant Nos. DE-FGO2-85ER-40214 and DE-FG05-91ER-40636 
and the Kansas Institute for Theoretical and Computational Science. 
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Figure Captions 

Fig. 1 Energy loss by emission of a parton in the presence of a source, marked 
by a circled "x" . 

Fig. 2 Space-time picture of the null plane uncertainty principle Ax~^Aq~ > 1 
The light-front time interval Ax~^ is bounded to be less than about 
V2La for ultra-relativistic processes. 

Fig. 3 Energy loss in a theory with a fundamental 4-point vertex; the momen- 
tum fractions are indicated by the lengths of the lines. The uncertainty 
principle bounds the formation rate of the fastest forming parton (for- 
mation rate k^/2xiP) but says nothing further about the loss occurring 
among the other two partons, provided the sum of their formation rates 
Fj = k^/2xiP is smaller than the first. 

Fig. 4a-d Ratio of the nuclear to Deuterium cross section for Drell-Yan contin- 
uum dileptons (DY) as a function of calculated by including the 
contribution only due to shadowing (dashed curve) and due to shad- 
owing and the maximum allowed value of energy loss (solid curve) for 
A=12, 40, 56 and 184, respectively. 

Fig. 5 The parameter a for bottomonium as a function of Xf calculated by 
including the contribution only due to shadowing (dashed curve) and 
due to shadowing and the maximum allowed value of energy loss (solid 
curve). 

Fig. 6a-d Ratio of the nuclear to Deuterium cross section for J/t/; production 
as a function of Xf calculated by including the contribution only due 
to shadowing (dashed curve) and due to shadowing and the maximum 
allowed value of energy loss (solid curve) for A— 12, 40, 56 and 184, 
respectively, using the transverse momentum distribution observed in 
the case of bottomonium. 

Fig. 7 Predictions for the limiting values of ratios of nuclear to Deuterium 
cross sections for T15 production as a function of xp. The data should 
lie above the curves. 
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Fig. 8 Limiting values of the ratio of Tungsten to Deuterium cross section for 
J/^l) production as a function of for different transverse momentum 
bins. The slopes of the curves at small transverse momentum are much 
smaller than the slopes at higher transverse momentum. Data in each 
transverse momentum bin should have a slope less than or equal to 
that of the plotted curve. The overall Xf independent normalization 
factor N is not predicted. 

Fig. 9 Ratio of Tungsten to Deuterium cross section for J/ip production as 
a function oi Xf including contributions of gluon fusion and quark- 
antiquark annihilation channels. The quark-antiquark annihilation con- 
tribution overtakes charmonium production for Xf > 0.65, producing 
the upturn in the curves. 
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